The Annals of Applied Probability 

2005, Vol. 15, No. 4, 2535-2552 

DOI: 10.1214/105051605000000575 

© Institute of Mathematical Statistics, 2005 

MONOTONE PROPERTIES OF RANDOM GEOMETRIC GRAPHS 
HAVE SHARP THRESHOLDS 

By Ashish Goel/ Sanatan Rai^ and Bhaskar Krishnamachari 

Stanford University, Stanford University and University of Southern 

California 

Random geometric graphs result from taking n uniformly dis- 
tributed points in the unit cube, [0, 1]'*, and connecting two points if 
their Euclidean distance is at most r, for some prescribed r. We show 
that monotone properties for this class of graphs have sharp thresh- 
olds by reducing the problem to bounding the bottleneck matching 
on two sets of n points distributed uniformly in [0,1]"*. We present 
upper bounds on the threshold width, and show that our bound is 
sharp for d= 1 and at most a sublogarithmic factor away for d > 2. 
Interestingly, the threshold width is much sharper for random geo- 
metric graphs than for Bernoulli random graphs. Further, a random 
geometric graph is shown to be a subgraph, with high probability, of 
another independently drawn random geometric graph with a slightly 
larger radius; this property is shown to have no analogue for Bernoulli 
random graphs. 

1. Introduction. Consider n points distributed uniformly and indepen- 
dently in the unit cube [0,1]'^. Given a fixed distance r > 0, connect two 
points if their Euclidean distance is at most r. Such graphs are called ran- 
dom geometric graphs, and are denoted by G^'^^ r), as in [24]. Classically, 
these graphs have been the subject of much study because of connections 
to percolation, statistical physics, hypothesis testing and cluster analysis. 
Further, random geometric graphs are better suited than more combinato- 
rial classes (such as Bernoulli random graphs) to model problems where the 
existence of an edge between two different nodes depends on their spatial 
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distance. As a result, random geometric graphs have received increased at- 
tention in recent years in the context of distributed wireless networks (such 
as sensor networks), see, for example, [12, 13, 14]; and layout problems as 
in [5, 16, 24]. Another area is cluster analysis, especially its applications in 
medicine, biology and ecology; these may be found in [10]. 

In applications such as distributed wireless networks, the connectivity 
of random geometric graphs is of interest. Gupta and Kumar showed that 
for d = 2, if 7rr(n)^ = (logn + Cn)/n, then as n | oo the graph is connected 
almost surely as n | oo if c„ j oo and is disconnected almost surely if Cn i 
— oo [12]. This result is remarkably similar to the corresponding result for 
Bernoulli random graphs (also known as Erdos-Renyi graphs). An instance 
of a Bernoulli random graph is obtained by taking n points and connecting 
any two with probability p, independently of all other pairs. This class of 
graphs is denoted by Gn,p- Erdos and Renyi [7, 8] showed that if p{n) = 
{logn + Cn)/n, then the graph is a.s. connected or disconnected as Cnl oo or 
Cn i — oo. For d = 2, Gupta and Kumar's result can also be obtained from 
Penrose's work on the longest edge of minimal spanning tree of G(A'„;r(n)) 
[23] . Connectivity of random geometric graphs for d= 1 was also studied by 
Godehardt and Jaworski [11]. Connectivity results under the /oo-norm may 
be found in [1]. 

In both random geometric graphs and Bernoulli random graphs, property 
thresholds are of great interest. To quote Bollobas [2]: 

One of the main aims of the theory of random graphs is to determine when a 
given property is likely to appear. 

Particularly interesting are thresholds for monotone properties, of which 
connectivity is a classic example. A seminal result of Friedgut and Kalai [9] 
states that all monotone graph properties have a sharp threshold in Bernoulli 
random graphs, and the threshold width is 6{e) = O {log /logn). They 
also demonstrated a monotone property with a threshold width of $7(1/ log^ n) 
and conjectured that this is tight [i.e., the best upper bound on threshold 
width is ©(l/log^n)]. Their upper bound on threshold width was improved 
to C'(l/log^~''' n) for all 7 > by Bourgain and Kalai [4]. 

The similarity of the connectivity threshold for random geometric graphs 
and Bernoulli random graphs led to the conjecture that all monotone prop- 
erties also have a sharp threshold in random geometric graphs (see [18, 20] 
for a more detailed discussion). For the d = 1 case, sharp thresholds for 
monotone properties are implicit in the recent work of McColm [20] , though 
he does not compute bounds on the width. The definition of sharp thresh- 
olds we use in this paper is the one used by Friedgut and Kalai and is based 
on the threshold width. The definition used by McColm is the one used in 
the text by Janson, Luczak and Rucihski [17], and is stronger than the one 
used by Friedgut and Kalai; we discuss this in more detail at the end of the 
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Introduction. The analysis of random geometric graphs is technically chal- 
lenging because of dependence of the edges. The triangle inequality implies 
that the event that points x and z are connected is not independent of the 
event {{x,y) and {y,z) are edges}. This is in stark contrast to the case of 
Bernoulli random graphs. Hence proof techniques that have been successful 
for Qn,p cannot be exploited in the case of random geometric graphs. 

Our results. We show that all monotone graph properties have a sharp 
threshold for random geometric graphs, thus resolving the above conjecture. 
In fact, the threshold width for random geometric graphs is much sharper 
than for Bernoulli random graphs. In order to state our results formally, we 
need to establish some notation and some definitions. 

We use the symbol to mean "distributed as," so that G~GW(A'„;r) 
means that G is picked from G^^\Xn;r). For ease of notation, we omit the 
superscript d in G^'^\Xn:,r) as the dimension will be clear from the context. 
The critical radius for connectivity is defined as Tc := (logn/irdn)^^'^ , where 
TTd is the volume of the unit sphere in W^. 

A graph property A is a set of undirected and unlabelled graphs. A prop- 
erty A is increasing if and only if 

GeA =^ {yG')[{V{G') = V{G) and E{G)<^E{G'))^G' e A]. 

Intuitively speaking, an increasing property is one that is preserved when 
edges are added to the graph. A graph property A is monotone if either 
A or A'^ is increasing. Without loss of generality, for the rest of the paper, 
we shall implicitly mean increasing properties when referring to monotone 
properties. 

If A is an increasing property, then for < e < 1/2, let r(n, e) = inf{r > : 
F{G{Xn',r) & A} > e}. Define further 5(n,e) := r{n, 1 — e) — r(n,e). A prop- 
erty is said to have a sharp threshold if 6{n, e) = o(l) for all < e < 1/2. 

Our main results are: 

Theorem 1.1. For every monotone property, the width 5(n,e) is 

for d = l. For d = 2, 

6{n, e) = 0{rc log^/^ n) = ©(log^/^ n/^), 
and for d>3, the width 

6{n,e) = 0{rc) = 0{\og^/'^ n/n^l'^). 
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Thus, all monotone properties have sharp thresholds. Observe that the 
width is much sharper than the threshold width for Gn,p- Moreover, we 
prove a stronger result: the graphs G{Xn]r) become subgraphs of G{Xn',p) 
for p > r, with hight probability (w.h.p.) as n | oo. Note only that this is 
not the case for Bernoulli random graphs — we shall make this precise in 
Section 2. 

For the lower bounds we have: 

Theorem 1.2. For d>2, there exists a monotone property with width 
6{n,£) = Q{{log£~^)^^'^n~^^'^). For d=l, there exists a monotone property 
with width 

Hence, we have a tight characterization of the threshold width for d=l, 
and our upper bounds are only a sublogarithmic factor away for d>2. 

The key idea is to relate the behavior of monotone properties to the 
weight of the "bottleneck" matching (to be defined later) of the bipartite 
graph whose vertex sets are obtained by distributing n points uniformly and 
independently in [0, l]*^. Sharp results on the "bottleneck" matching weight 
are implicit in the work of Leighton and Shor [19] for d = 2 and Shor and 
Yukich [26] for d>3. We repeat them here for convenience. 

Theorem 1.3 ([19, 26]). Consider the bipartite graph on 2n points, 
where each set of n points is distributed uniformly and independently in 
the unit cube [0,1]"^, and independently of each other. If M is the length of 
the bottleneck matching, then w.h.p. as n] oo: 

rG(relogi/^n), ifd = 2, 
le(re), ifd>3. 

In Section 4 we present our own proof of the bound for d>3. The proof 
for the d>3 case in Shor and Yukich [26] invokes results from polyhedral 
geometry. We present a simpler proof that relies only on the properties of 
order statistics and Chernoff bounds. We prove that for d = l, the bottleneck 
matching is 0{y^log£^^ / y/n) with probability 1 — e. Moreover, our results 
also hold for any ^p-norm when p > 1, and not just under the Euclidean 
norm as in the setting of this paper. We omit the details, as they require 
straightforward modifications of the proofs given herein. 

It might seem curious that we do not report the dependence on e in some 
of our bounds on the threshold width. This is because the results of Shor and 
Yukich as well as those of Leighton and Shor are high probability results: the 
bottleneck matching length is 0(rclog^/^n) in two dimensions and 0(rc) in 
higher dimensions not just in expectation but with probability 1 — o(n~^) 
for some (3 > 0. Hence, in asymptotic notation, our upper bound on 6{n,e) 
in two and higher dimensions does not depend on e as long as e = r2(?i~^). 
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Related work. There is a vast body of literature that is directly related 
to this paper. It would require a survey paper to even mention the salient 
results with any degree of honesty. We can only point the reader to the book 
by Penrose [24], the papers by Gupta and Kumar [12, 13, 14] and the paper 
by Shakkottai, Srikant and Shroff [25]. We note here that our techniques 
imply a sharp threshold for the coverage problem as discussed in [25] , which 
is not a graph problem. We omit the details. The theory of Bernoulli ran- 
dom graphs is covered in the books by Bollobas [2] and Janson, Luzak and 
Rucinski [17]. For some results on matchings in a similar context, see the 
paper by Holroyd and Peres [15] and for some results on covering algorithms 
see [3]. Sharp thresholds for random geometric graphs were conjectured in 
[18] and [20]. Muthukrishnan and Pandurangan [21] obtained asymptotically 
tight thresholds for connectivity, covering and routing-stretch in d dimen- 
sions using a new technique called bin-covering. 

Additive versus multiplicative thresholds. In this paper we are primarily 
concerned with bounding the threshold width of a property, along the lines of 
Friedgut and Kalai [9]. Informally, this corresponds to proving sharp "addi- 
tive" thresholds. As we mentioned earlier, the notion of sharp thresholds pre- 
sented in [17] or in [20] is stronger in that they require 5(?l,e)/r^(?^) = o(l). 
Informally, this corresponds to "multiplicative" thresholds. We observe that 
our Theorem 1.1 also yields sharp thresholds in this stronger sense, provided 
the threshold radius is high enough. More precisely, if 11 is a monotone prop- 
erty and ^^(^^) is its threshold radius, such that: 

when d>3, 
when d = 2, 
when d=l, 

then n also has a sharp threshold in the sense of Janson, Luczak and 
Ruciiiski [17]. 

Plan of this paper. We first establish the relationship between monotone 
properties and bottleneck matchings and prove the upper bound in Section 2. 
In Section 3 we furnish the lower bounds. In Section 4 we discuss the upper 
bound for d>3, and in Section 5 for d= 1. We conclude in Section 6 with 
some open problems. 

2. Bottleneck matchings and monotone properties. Recall that in a bi- 
partite graph with vertex sets Vi and V2, a perfect matching is a bijection 
(/) : Vi — > V2, such that each v is adjacent to (p{v). Thus a perfect matching is a 
disjoint collection of edges that covers every vertex. If the graph is weighted, 
then we define the weight of the matching as the maximum weight of any 



Tc = o(rn(n)) 
rc = o(rn(n)/log^/^n) 
y/nrii(n) — > 00 
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edge in the matching. A bottleneck matching is the perfect matching with 
the minimum weight. 

Let Si and S2 denote two sets of n points each, where the points are i.i.d., 
chosen uniformly at random from the set [0, l]'^. Form the complete bipartite 
graph on (Si, 5*2 ) and let the weight of an edge be the Euclidean distance 
between its endpoints. Let Mn*^^ denote the bottleneck matching weight of 
this graph. We omit the dimension d where it is clear from the context. 

We first link the weight of the bottleneck matching with a containment 
property on random geometric graphs. We shall write G C G' to mean that 
the graph G is contained in the graph G', that is, is isomorphic to a subgraph 
of G'. 

Lemma 2.1. Suppose P{M„ > 7(n)} <p for some function 7(n) and 
some constant p. For any radius r, consider independent random samples 
G~G(A'„;r) and G' ^ G{Xn]r + 2'y{n)) in d dimensions. Then, P{G C 
G'}>l-p. 

Proof. Let V represent the set of points in graph G and V' the set of 
points in graph G' . Let (j) denote the bottleneck matching between V and 
V'] then Mn is the weight of this matching. Suppose (it,f) G E{G), that is, 
11^^ — II2 ^ Then, using triangle inequality, 

\\(j){u) - (l){v)\\2 < ||0(n) - u\\2 + \\u - v\\2 + \\v - 4>{v)\\2. 

But ||</'(u) — ti||2 and ||^(i') — v\\2 are both at most M„, and hence ||</'(ti) — 
0(^^)||2 ^ 2M„ + r. If M„ < 7(n), then the mapping (/> establishes that G C G' , 
and hence F{GcG'}>l-p. □ 

The main result linking monotone properties to bottleneck matchings is: 

Theorem 2.2. //P{M„ > 7(n)} <p, then the ^/p-width of any mono- 
tone property in d dimensions is at most 2'y{n). 

Proof. Let p = e^, so that P{M„ > 7(n)} < e^. Let n be an arbi- 
trary increasing monotone property. Let = r(n,e), ru = r^ + 27(71). Let 
G ~ G{Xn;rL), and G' ~ G(A'„;rc/), and define q := F{G' ^ U}. Since G is 
independent of G' , ¥{G £ll,G' ^ H} = e ■ q. The monotonicity of n implies 
that if G G n and G' ^ H, then G G' . This means that P{G (^G'}>e-q. 
By Lemma 2.1 above, P{G G'} < p, so that we must have q <£. But then 
r(n, 1 — e) < r(n, 1 — q) = ru, so that 6{n, e) <ru — r^ = 27(72). □ 

With Theorem 2.2, the upper bound theorem follows with very little more 
work: 
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Proof of Theorem 1.1. Leighton and Shor [19] show that 

with probabihty at least 1 — n^'^, for some k> 0, and Shor and Yukich [26] 
show that 

71/4"^)= e(rc) ford>3, 

with probability at least 1 — n~'^' , for some k' > 0. Hence Theorem 2.2 implies 
that 6{n,e) = 0{rclog^^^ n) for d = 2 and 5{n,e) = 0{rc) for d>3, for any 
constant e > 0. In fact, the bound on 5{n, e) holds for any e = Q{n~^), where 
c > is a constant. 

In Proposition 5.1 (see Section 5), we show that for d = l, 

F|M«<-^|>l-exp(-c/32), 
for some c > 0. By applying Theorem 2.2 with e = exp(— c/?^) we obtain 



5{n,£) = 0\ 



llogt 



n / □ 



Remark 1. We have in fact shown that G{Xn;r) is a subset of G{Xn] f') 
w.h.p., when r' = r + o(l). The corresponding result does not hold for 
Bernoulli random graphs. To see this suppose that G and G' 

For G cG' , every edge in G must exist in G' . Hence, for M = C^), q = 1 —p 
and Q = 1 — P, and a given matching (p: 

no C G' under </>} = E (^) P^'q^'-'' 
K=0 ^ ^ 

^ ' M - K\ ryK+LriM~K~L 



L=0 

M 



pK^M-K 



K=0 

M-K 



E(^)(p^)V-^(^'+g)^^-^ 

{pP + q) 
{p{l-Q)+q) 
{l-pQf'. 



K=0 
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Choose p = 1/4, and P = 3/4. As there are n! matchings: 



•{GCG'}<„!exp(-!*i^). 




The last expression goes to zero as n | oo. Hence, even in this extreme case 
when P — p = l/2, we do not have Gn,p C Gn.p with constant probability. 

3. The lower bounds. We now present examples of monotone properties 
to show that our bounds are tight in the d= 1 case and within a subloga- 
rithmic factor for d>2. 

For the d = l case, we consider the property 11, defined by 

\V\ 

GgIT <^=^ mindeg(n) > — — , 
u&v 4 

where V is the vertex set of G. Let xi, . . . , x„ be the n uniformly distributed 
points in [0, 1]; these are the vertices of the graph G. Let denote the ith 
order statistic. We have the following two estimates: 

Lemma 3.1. IfO<e< 0.5e~^^/^, then for property U: 
r{n, 1 - 

Proof. Let u= 1/4 + A, where A > is to be determined later. Pick 
G^G{Xn',u), and let the vertices be xi,...,Xn- Then dis- 
tributed uniformly in [0, 1] . Observe that 

P{x(,/4) > n} > e =^ p|deg(x(i))<^|>e 

=^ F{G ^ n} > e, 
where in the first implication we have used the fact that deg(a;(i)) < n/4 44> 

X{n/4+i) - > u, and that X(„/4+i) - = X(„/4). 

Now, P{x(„/4) > u} = P{Bin(n, u) < n/4}, and for some suitably large no, 

P{Norm(0,l) <-V«^A}>2e =^ F{Bin(n, n) < n/4} > e, 

whenever n > no, by the Normal approximation to the Binomial. 

Put A = /3/V^ for (3 = V61og(0.5/e)/ll. Then for < e < 0.5e-^^/^ we 
have P >2. With a little bit of work, one can see that (3'^/ 2 < -4/3/3 - 
log2e. Since x > logx for x >1, the last inequality implies that /?^/2 < 
log(3/3~-^/4) — log(2e). Observe that any /? > 2 satisfies 

4/? - /3 /33 
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so that 

/32 / 1 1 \ 1 



2 ^^°^U"^J+'°'2e' 



or 

32 N 



(/3-'-/3-')exp(^-^)>2e. 



But then by Theorem 1.4 of [6], we can conclude P{Norm(0, 1) > /3} > 2e. 
This shows that 



and since 11 is increasing, this means that 



Lemma 3.2. For property U: 

/ 1\ 1 c 

r n, — < — I —, 

\ '2/ - 4^ V^' 

for some fixed constant c > 0. 

Proof. Suppose G ~ G(A'„, /). For any u G ^(G), write degL(n) for the 
number of points to the left of u and adjacent to it, and similarly let degj^(M) 
stand for the number of right neighbors. Note that if deg(a;(i)) > n/4, then 
necessarily deg(a;(j)) > n/4 for all 1 < i < n/4. With this observation we 
have: 

P{G^n}=p| U |deg(x(,))<^|| 

< p||deg(x(i)) < ^ or deg(x(„)) < ^|| 
U |degL(a;(.)) < ^1 

.n/4<i<3n/4 ^ ^ ) 



+ p| U |degR(x(,))<^U 
< 2p|deg(x(i)) < ^1 +np|deg(x(,)) < 



n 



(1) (2) 
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To bound the first term (1), first observe that by arguing as in the 
last lemma P{deg(x(i)) < n/4} = P{Bin(n, /) < n/4}. By applying Chernoff 

bounds, we can find a Ci > so that P{Bin(n,/) < n/4} < e"*^?/^, when 
/ = l/4 + Ci/V^. 

To bound the second term (2), again apply Chernoff 's bounds to find 
C2 > 0, such that for / = 1/4 + C2I ^f^, 



np|deg(a;(i)) < '^j = "-Ip|Bin(n, < ^| < nexp 



n 
'32. 

so that for n large enough this term is overwhelmingly small. Therefore, for 
c > max(Ci , C2) , and / = 1/4 + c/ y^, we have 

P{G^n}<e-^'/2 + nexp(^-^), 

so that r(n, 1/2) < 1/4 + c/-y/n, for a suitably chosen c. □ 

Lemmas 3.1 and 3.2 show that for the graph property 11 defined by 

\V\ 



GgH <^=^ mindeg(u) > , 
u£V 4 



we have 



rn(n, 1 - e) > 7 + ^hen < e < 0.5e"^^/^ 

4 V 2n 



rn(n,-) <- + ^. 



Hence we have shown the d > 2 case of Theorem 1.2: 

Proof of Theorem 1.2, lower bound for d=\. Immediate from 
the last two lemmas. □ 

Theorem 3.3. For d > 2, there exists an increasing property H such 
that for < e < 1/2, the threshold width satisfies 

Proof. Let n be the property that the graph is complete. This is triv- 
ially a monotone property. 

Suppose < e < 1/2. Let u := y/d{l - 2A+) [see Figure 1(a)], for A+ such 
that 

< A+ < -^[min(^,log2/2)]^/'', 

and pick G ~ G{Xn;u). Fix a pair of diagonally opposite corner cubes with 
side A_|_, and consider the event that there is exactly one point in each. If 
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this happens, then the graph is not complete, since the points are more than 
u apart. Hence: 

P{G ^ n} > 2 ) (^+)^ • 2 • (1 - 2A^)"-2, 
since A+ < (log 2/(2?i))^/'^. Thus, for large enough n we have 

which implies that 

P{G'^n}>(^)(Ai)2>e. 

The last inequality follows because we chose A'^ < 7i^^ x min(\/2e, log 2/2). 
Therefore, 

(1) r{n,l- e)>u>Vd{l 



n 



l/d 



Now we shall bound r{n) above. To this end set / = — 1 + (1 — 4A_)2 
[see Figure 1(b)]. Now suppose that A_ = (loge~"'^)^/'^/(4n"'^/'^), and pick 
G^G{Xn]l)- Using elementary geometry it is easy to see that, if none of 
the n points lies in any of the 2*^ cubes of side 2A_ at the corners of [0, l]"^, 
then the graph is complete. Hence, for n large: 

P{G G H} > (1 - 2'='(2A„)'^)" > exp(-n(4A„)'^), 

where we have used the fact that l — x> , when x > 0, so that (1 — x)" > 
e~"^, if X < 1. Since exp(— n(4A_)'^) = e, by our choice of A_, it must be 
that 



(2) ^n(^)<^ = ^^-l+(l- ^'°',v7^' ) 



Putting (1) and (2) together: 
6{n, e) = r„(l - e) - r„(e) > M - Z 



n 



Vd 



ce^/d\ / 2/loge-iy/'^ 1 /loge-i\ ^/"^ 
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(a) 



(b) 



Fig. 1. Definition of A+ and A- 



Vd 



Vd 



n 



"^2dV n ) V ~2V n 

1 4 / logg-^ y/'^/ 1 / logg-^ y/^y 

logV'ie-i 



n 



□ 



Observe that for any k > constant, \ie = n~'^, our lower bound for d>3 
matches our upper bound on the threshold width, and is only a factor of 
0(log^''^n) away for d = 2. For any constant e, the difference between the 



bounds is ©(log 



l/d 



n) 



and 0(log^/^n), respectively. 



4. Bounding the bottleneck matching for d > 3. We now present a sim- 
pler proof of Mn = Q{rc) when d>3 than the proof presented in [26]. We 
emphasize that even though we also recursively subdivide the cube, our 
principle is different. In our proof, at the end of each step, the points are 
distributed uniformly in each subcuboid. This requires careful choice of the 
cutting plane and a linear transformation based on order statistics. This, 
however, permits us to bound the matching length via Chernoff bounds, as 
opposed to estimating the aspect ratios of rectangular solids as in [26]. 
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The basic idea is to divide the unit square into n equal boxes. Given n 
points distributed uniformly on the square, move the points so that there 
is roughly a single point in each box. Now consider the two samples of 
red and blue points. Apply this process to both samples. Let A be the 
maximum distance by which a red point is shifted, and similarly, let A' be 
the maximum shift for any blue point, along any coordinate direction. Then 
the triangle inequality tells us that the bottleneck matching is less than 
V5(A + A' + n-i/"^). 

We shall now use this idea to bound the bottleneck matching in [0, 1]"^. 

To move each point into its unique box we follow a recursive process. 
We shall provide only an informal description here, relegating the details to 
the Appendix. Moreover, for simplicity, we shall suppose n to be a power 
of 2. First divide the square by drawing a vertical line so that there are 
exactly n/2 points in each half. Transform the x-coordinates of the points 
in each half so that they are uniformly distributed in [0,1/2] and [1/2,1]. 
Now repeat the process along the y-axis for each rectangle, and then along 
the z-axis and so on. Repeat when all coordinate axes have been done once, 
and so on. One can carry this process for log n steps so that there is exactly 
one point in each box. However, for d > 3 it is better to stop at the jth step, 
where j < logn. Choose j = log2(n/ logn)] . Then the side of the box 
and hence the shift thereafter is < . With this observation we can now 
show 

Proposition 4.1. If Mn is the weight of the bottleneck matching on a 
geometric random bipartite graph on 2n points in [0, 1]"^, for d>3, then for 
any j3> 1, we can find a constant Q > such that 

P{M„>Cd/3rJ<-^, 

so that Mn = 0{rc) w.h.p. 

Proof. To estimate M„, we compute the total shift experienced by an 
arbitrary point. To this end, we shall find the shift along each axis, and 
so shall concentrate on one coordinate at a time. Let xi, . . . ,2;^ denote the 
coordinates. We regard a step as one cycle in which divisions along all axes 
have been completed, according to the scheme described above. Therefore, 
if a step begins with a d-dimensional cube of side / containing n points, by 
the end of the step, the cube has been divided into 2'^ cubes of side //2, with 
n/2'^ points each. 

Let Ui = n/2'^(*~^^ denote the number of points in a subcube at the be- 
ginning of the ith step, and let li = 2~*+^ denote the length of the side of 
the cube. Lemma A. 2 implies that for any point in the left half of such a 
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subcube, the shift 5) in the x^-direction experienced during the ith step 
satisfies 



Therefore, if 6i is the total shift suffered by a point during the ith step 



Now fix a /3 > 1, and choose such that ^yf ■n/2^^~^^^'^~'^^ = 0^ logn. Observe 
that with this choice, 7^ is decreasing with i only when d> 2. Let A be the 
maximum total shift experienced by any point. Then it must be that 



which follows from taking the union bound over all n points, and all d 
coordinates, and the fact that after j steps, we divided a given coordinate 
at most j times. Notice that after j = logn steps, the side of the subcube 
reduces to 2"^°^""'"^ = 0(l/n), and therefore subsequent shifts cannot move 
the point by more than 0{l/n). Hence we can halt the subdivisions after 
logn steps, knowing that the matched point is already within 0(l/n). Hence, 




'i 



for any 7i > 0. 






L 1=1 J 

However, with a little bit of work, one can see that 







Recall that {log n/n)^/'^ 




so that we have just shown 
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Therefore, we have 

n|A| > 2Vdf3.'fr.}<Fl\A\ > Vd±j] < 



i=l 



After j steps the side of the cube is 2 and hence if we arbitrarily move 
points within the subcube, the extra shift is at most Therefore, if we 

choose Cd to be any constant larger than \/^ + 2\/d7ry'^, we get |A| < Cdl3rc 
with probability at least 1 — n^~f^^ . □ 

We note in passing that for the above method to provide a bound in the 
d = 2 case, one must proceed for logn steps, so that there is only a single 
point in each box. However, in this case, one only gets an 0{rclogn) bound, 
which is off by log^^^n from the sharp bound in [19]. 

5. The bound for d= 1. Given n points uniformly distributed in [0, 1], 
follow the recursive division procedure described in the last section. In this 
case, at each step the number of points decreases by half. Therefore, we 
obtain a stronger result: 

Proposition 5.1. For any /? > 

P{M^,^) > (3/^} = 0(exp(-c/32)) 
for some positive constant c. 

Proof. For the sake of simplicity we assume that n = 2^ for some /c G N; 
this makes no difference to the proof except for simplifying some expressions. 
In the ith step there are 2* sets of n/2* points each. Therefore, if 6i is the 
shift of an arbitrary set, then for Pi > 0, by Lemma A. 2: 

p{|<5.|>^}<2exp(-2*/?2), 

so that the maximum shift A of any point satisfies 

p{max|A| > ^1 < 2^|exp(-2^/3f). 

Choose the /3j's such that 2*/?? = /?q + i, for some Pq. Then J2i A < for 
some suitable constant K > 0. Taking P = KPq, we get 

P 



p|M«>-^|<c'exp(-c/32), 



for some constants c, c' > 0. □ 
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6. Conclusion. We have shown that all monotone graph properties have 
a sharp threshold for random geometric graphs. Moreover, this threshold is 
sharper than the one for Bernoulli random graphs. 

We have a sharp result for d = 1. For the d > 3 we believe the upper 
bound of 0{rc) to be actually tight. For the d = 2 case we believe the upper 
bound to be 0{rc) as well, though this cannot be obtained via bottleneck 
matchings. 

APPENDIX: ESTIMATING THE SHIFT IN EACH RECURSIVE STEP 

To establish a bound on the amount by which each point is moved, we 
must examine the "shrinking" and "stretching" process formally. For sim- 
plicity we concentrate on the d = 2 case. Ignore the y-coordinates. Then we 
have n i.i.d. Unif [0, 1] points xi,. . . , Xn- It is well known [22] that the k small- 
est points are i.i.d. uniform in [0, x*^*^"*"^^), and that the n — k largest points 
are i.i.d. uniform in (x^'^^l], where x^*^ is the ith order statistic. Suppose 
that n is even — the analysis below applies mutatis mutandis when n is odd. 
Set 5i = x("'/2+^) — 1/2 and (5,- = 1/2 — x^"/^). Then transform the points as 
follows: 

x^*) — }-^-^—r, forz = l, ^ 



X 



l/2 + <5;' '■■■'2' 
-(l-x«)— for i = - + l,...,n. 



1/2 + 5r 2 

This transform leaves the smallest n/2 points uniformly distributed in [0, 1/2] 
and the largest n/2 points uniformly distributed in [1/2, 1]. This process is 
now repeated [logn] times alternating the x- and y-coordinates. The max- 
imum shift at any step is not more than \5i\ for the points on the left and 
not more than \5r\ for points on the right. We shall use 5 = max{5;,(5r}. 

Let Xt be the number of points in [0,t] prior to the transformation. The 
following result is immediate: 



Lemma A. 1. For 0<7< 1/2, 

P{|5|>7}<2p|xi/2+-,<^}. 



To bound the last probability, observe that Xt is just the sum of n i.i.d. 
Bernoulli's that are 1 with probability t. Hereafter /J > is some constant. 

Lemma A. 2. The shift 5 of any point in the recursion step satisfies 

P{|5|>7} = O(exp(-n'(7/0')), 

where n' is the number of points in the subcuboid being divided, and I is 
length of the side that is being divided. 
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Proof. The proof is a straightforward apphcation of Chernoff 's bound. 
Assume wlog that / = 1: 

P{|5|>7}<2p|xi/2+,<^} 

= 2p|xi/2+^ < Q + 7) - n'7| 

< 2exp(-nV) since 7 < 1/2. □ 
Generahzation to the d>3 case is straightforward. 
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